Vortex Charging Effect in a Chiral p x ± ip y -Wacve Superconductor 
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Quasiparticle states around a single vortex in a p x ± ip^-wave superconductor are studied on 
the basis of the Bogoliubov-de Gennes (BdG) theory, where both charge and current screenings 
are taken into account. Due to the violation of time reversal symmetry, there are two types of 
vortices which are distinguished by their winding orientations relative to the angular momentum of 
the chiral Cooper pair. The BdG solution shows that the charges of the two types of vortices are 
quite different, reflecting the rotating Cooper pair of the p x ± ip^-wave paring state. 



I. INTRODUCTION 



a 
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The discovery of many types of superconductors from heavy fermion compounds to high-T c cuprates has driven us 
5h ■ to study a large variety of new physics beyond the standard BCS theory for conventional s-wave superconductors. 
The study of the unconventional superconductivity was stimulated by the discovery of superfluid 3 He in which a 
spin triplet p-wave state realizes. Unlike the conventional s-wave state, the p-wave state has both spin and orbital 
degrees of freedom. This is the most pronounced feature of the unconventional superconductors observed in their 
thermodynamics and impurity effects or detected by tunneling spectroscopy, NMR, /iSR measurementSpand so on. 

Sr2RuC>4 is the first layered perovskite compound showing superconductivity without CuC>2 planes .a The recent 
experimental and theoretical studies have indicated that the superconducting pairing symmetry of S^RuC^ is not a. 
simple s-wave. The absence of a Hebel-Slichter peak in NQR0 and the sensitivity of T c an non-magnetic impuritiestl 
point towards unconventional pairing. The indication of broken time reversal symmetryp observed in /xSR measure- 
ments, gives a strong argument for the unconventional parine. state. The Knight shift experiment shows that the 
spin susceptibility is not affected by the superconducting state,Q which is the strong evidence of a spin-triplet pairing. 
Sigrist et al. suggested that & p x ± ip y -w&ve state^-which breaks the time reversal symmetry in a tetragonal crys 
field, is the most likely pairing state for S^RuC^.B The line node behavior is reported in the latest experiments 
related to the low temperature thermodynamical peasurements, such as specific heat and NMR T x . An orbital 



dependent superconductivity!!!] and gap anisotropyudH were suggested to understand the line node behavior. Here 
we focus on the p x ± ip^-wave pairing state, since this representation is the simplest and essential form. We will see 
t— I ■ rich physics of this chiral state. 

The most intriguing character of the p x ± ip y -wave state is that the Cooer pair has ±1 angular momentum, i.e. the 
pair electrons are rotating. This property is similar to that of the A phase of the superfluid 3 He. Due to the violation 
of the time reversal symmetry, we have two types of vortices, one of which is in the same direction to the angular 
momentum of the rotating Cooper pair, and the other is in the opposite direction. The rotating pair shows up in 
quasiparticle states around a vortex core. _ 
In this letter we present an interesting new physics related to a vortex. We focus on a vortex charging effect .Ej 



■ It was pointed out that for_an s-wave superconductor the vortex charge is proportional to the slope of the density 



of states at the Fermi level 



14 



Hayashi et al. has proposed that the vortex charge is determined by the quasiparticle 



structure in the vortex corel!3 rather than the slope of the density of states. Very recently, it has been reported that 



Chern-Simons terms lead to a fractional vortex charge for the p x ± ip y -w&ve stateJia Thus the origin of the vortex 
charge is still controversial. To make this point clearer, we investigate the vortex charging effect in the chiral p x ± ip y - 
wave state, concentrating our attention on a microscopic origin of the vortex charge. This is demonstrated by solving 
the BdG equation self-consistently, including both charge and current screenings. This is the first full self-consistent 



BdG study of a single vortex. We would like to show how the rotating Cooper pair shows up in the vortex charging 
effect. 

This paper is organized as follows. In Sec. II, we present our formulation for calculating the quasiparticle states 
near the vortex core based on the BdG equation. In Sec. Ill, we show numerical results. We give summary and 
discussions in Sec. IV. 
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II. FORMULATION 



Let us begin with the following BdG and gap equations for the p x ± ip y -wave state 
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The detailed derivation is given in the Appendix. Throughout this paper, we use H = 1 and k^ = 1 units, e > is 
the electron charge. E n is the energy eigenvalue for the superconducting quasiparticle, and D± = d x ± id y is used. 

and A are the scalar and vector potentials, respectively. V p > and f2 are the p-wave attractive interaction and 
volume of the two-dimensional system, respectively. E c is a cutoff energy, and /i is a chemical potential. A± is the 
order parameter for the p x ± ip y -wave. For simplicity, we assume that the superconductor is basically two-dimensional 
and has a cylindrical Fermi surface. We note that the BdG solution has the following time reversal relation: 



{u-E n ,v~E n } = {v* En ,u* E J. 
The solution of the BdG equation determines the two-dimensional electron and current densities: 



n(r) = 2 
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\u n (r)\ 2 f(E n ), 
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By using the relation (g), both densities have been expressed in terms of the u n amplitude. f(E n ) is the Fermi 
distribution function and the E n summation in Eq. (|3|) runs both negative and positive regions. The scalar and 
vector potentials obey the Maxwell equations: 



V 2 A (r) = -^[n -n(r)] 
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We have introduced a layer spacing d to convert the area densities into the volume ones. To satisfy charge neutrality 
we introduce a uniform ng as the density of positive background charge. We have taken the origin of the coordinate 
at the vortex center. 

In the bulk region p x + ip y and p x — ip^-wave states are degenerate. In this letter we will choose one of the two 
degenerate states, the p x + ip y -w&ve state as a dominant component. However, the other component (p x — ip y ) is 
admixed with the bulk state (p x + ip y ) close to the vortex core.E3 Therefore our formulation includes both p x ± ip y 
components. The p x + ip y -w&ve state has a +1 Cooper pair phase winding, since □+ in Eq. (Ill) is expressed as 
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(5) 



Here r and (j) are the two-dimensional polar coordinates. A single vortex produces an additional phase winding ±1. 
Therefore there are two kind of vortices: the vortex winding direction is parallel or anti-parallel relative to the Cooper 
pair winding. As we will see later, the former vortex is charged up, while the latter is not. We name the former and the 
latter C(charged)-vortex and U(uncharged)-vortex, respectively. The C-vortex has two fold winding +2 effectively. 
On the other hand, the phase winding is zero for the U-vortex, since the winding of the Cooper pair cancels the 
vortex one. Due to the rotational symmetry of the system, angular momentum is a good quantum number. The BdG 
solutions are then classified by the angular momentum. To keep the rotational symmetry, there is only the following 
combination of the order parameters for the C and U- vortices: 
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(C- vortex) A+(r,0) = A+(r)e** A_(r,0) = A_(r)e l3 *, 
(U-vortex) A+(r, 0) = A + (r)e"^, A_(r, 0) = A_(r)e^, 



(6a) 
(6b) 



where A + and A_ are the dominant and admixed components, respectively. Hence the wavefunction u„ couples with 
v n in the following angular momentum spaces:E3 



(C — vortex) u l n «-> v: 
(U — vortex) u l n 
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where the superscripts represent the angular momenta of u n and v n . The GL calculation showed that the U-vortex 
is energetically favored.Eil However, in a real sample we can expect two types (p x ± ip y ) of domains, so that there 
are both C and U-vortices in the presence of an external magnetic field. We solve the two-dimensional single vortex 
problem on a disk of radius R. Since we treat a cylindrical system, it is convenient to expand wave functions by the 
following base functions: 



u n (r)\ = 
v n{r) J 
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Here Ji is the l-th Bessel function, Zij is the j-th zero of J/ and R is the radius of the system. I' in Eq. (|9|a) takes 
1 — 2 (I) for the C- vortex (U-vortex). 

In the practical numerical calculation, we use no-dimensional quantities. We express them with bar as follows: 
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Here Ao (unit of energy) is the magnitude of the order parameter at T = 0. £,o — vf/Aq (unit of length) is the 
superconducting coherence length. Using these non-dimensional quantities, we obtain the following set of equations: 
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V 2 a (f) = --— [no-n(r)], (10g) 

V 2 d(f)=-ff) 2 ij(f). (lOh) 

Here Atf and Al are the Thomas-Fermi screening lenth and London penetration depth, respectively. They are given 
by 

A ^ F = 4vre 2 7V(0) = 4^' 
2 c 2 md 
47T e^no 

where iV(0) is the density of states at Fermi energy. For two-dimensional layer system, it is given by iV(0) = mjixd. g 
in Eq. (|10p) is a non-dimensional coupling for the p-wave defined as g = V P S1, / (A^t; 2 ) . We have three no-dimensional 
parameters: /cf£o, Atf/£oj and Al/Co- 

In our numerical calculation, we obtain A±(f), ao(f), and a(f) by solving Eq. (|l^) iteratively within a Gygi and 
Schliiter methodB3 In the self-consistent calculation we fix the total number of electrons to the normal state value by 
adjusting ft. 



III. NUMERICAL RESULTS 



Let us discuss our BdG self-consistent solutions. First we show the order parameters in Fig. 0(a). The admixed 



3 



component A_ is induced around the vortex core as expected. It shows_t.hc asymptotic behavior A_ oc r (A_ cx r 
) for the U- vortex (C- vortex), which is consistent with the GL result E3 |A_| for the U- vortex is larger than the 
C- vortex one. This indicates that the U-vortex gains much condensation energy. In Fig. |(b|Lwe show several energy 
eigenvalues of the BdG equation. We notice that both vortices have zero energy bound states.E3 The appearance of the 
zero energy states is the consequence of the symmetric property of the BdG,-equation [Eq. (j^)]. As we discuss later, 
the bound state for l u = is very important to the vortex charging effect JISEj where l u is the angular momentum of 
the wavefunction u n . Be and Bu in Fig. l](b) represent the bound state. Next we show the charge density around 
the vortex core. As in Fig. ||(a) large charge density appears in the vortex core.12.3 The induced electric field is 
screened as we go far from the vortex center. Figure ^(b) is the spatial dependence of the electron density at various 
temperatures. At T = the electron density is suddenly decreased in the core region, which results in the vortex 
charge. With the increase of temperature, the electron density becomes uniform and the charge density is reduced 
accordingly. Contrary to the C-vortex, the electron density is almost uniform at all temperatures for the U-vortex, 
so that the vortex charge is very small in this case. 

Let us explain the microscopic origin of the vortex charging effect. First we discuss the C-vortex case. There are 
two contributions to the electron density n. One is from the bound states, and the other is from the extended states. 
We name the former n^, and the latter tie [n — + ng). In general only the l u = states contribute to the local 
electron density at the vortex center.. In the microscopic study we find out the important role of the bound state with 
zero angular momentum (Z„ = 0).li3o At zero temperature only the E n < states are effective in Eq. (^). Therefore 
the bound state Be in Fig. |l|(b) cannot contribute to the electron density, since this state is unoccupied [Fig. ||(a)]. 
«b is then suddenly decreased in the core region as in Fig. ||(b), and the total electron density is decreased close to the 
vortex center, resulting in the finite vortex charge. At finite temperatures the contribution from the Be state comes 
out [Fig. ||(a)] due to the finite Fermi distribution function in Eq. (||) for E n > 0. Correspondingly, tib increases in 
the core region as shown in Fig. ||(b). At T = 0.3Ao, which is larger than the energy of Be, the electron density 
is almost uniform due to the contribution from the Be state. The vortex charge is then reduced strongly at high 
temperatures as shown in Fig. |^(b). 

Next we discuss the U-vortex. Contrary to the C-vortex, the Bu bound state in Fig. [|(b) can contribute to the 
electron density at all temperatures, since this state is located at zero energy. As the consequence, the electron density 
for the U-vortex is almost uniform. In the U-vortex case we find very small temperature dependence in the electron 
density, since the Bu state is located at zero energy and participates in the energy excitations with the same rate 
at all temperatures. We note that the electron density for the U-vortex is very similar to that for the C-vortex at 
T = 0.3Ao, where the Be state is occupied. We conclude that the appearance of the vortex charge depends on the 
position of the l u = bound state relative to the temperature. In the conventional s-wave case, the vortex charge 
always appears at sufficiently low temperatures,E3 which is similar to the C-vortex result. 
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Let us mention the effect of charge screening. We show the total electron density n and contribution from extended 
states He for two cases, where charge screening is taken into account or not (see Fig. In the no charge screening 
case, the contribution from extended states tie is small at the vortex core. The total electron density n is then 
decreased, resulting in larger vortex charge. Note that the scalar potential in Eq. (|Io|g) is zero (oq = 0) for no charge 
screening (At — * oo). Therefore the scalar potential does not work to screen the vortex charge. On the other hand, 
when the charge screening turns on, quasiparticles respond to the vortex charge via the scalar potential do in the BdG 
equation (JToja.)- (|To|c) . tie is then increased to cancel the vortex charge. In contract to the extended state, electron 
density of the bound states iib does not exhibit difference between the charge screening and no screening cases. Since 
the bound states are localized around the vortex core, their wave function hardly modulate, while the extended states 
do easily. Hence the extended states quickly respond to the electric field and they screen the vortex charge. However, 
there still remains the substantial vortex charge. I— ■ I— ■ 

We discuss the effect of magnetic field. The Chern-Simons physicslla and Bernoulli effected take place in the presence 
of magnetic field. To see this effect we compare the two results, where a is present or ignored (no magnetic field). The 
scalar potential do is taken into account for both cases. Our result shows that the effect of the magnetic field is too 
small to be seen. Thus we find that the vortex charge is mainly determined by the microscopic quasiparticles structure, 
which reflects the phase winding of the chiral pairing. Figure ^| shows a schematic picture of the energy spectrum 
of vortex bound states for various winding pairings. The lowest-lying state appears at the symmetric position of l u 
and l v due to the relation (0), where and l v are the angular momenta of u n and v n , respectively. For the C- vortex 
the bound state energy is positive for l u — 0, while it is negative or zero for the U-vortex. Therefore the C-vortex is 
always charged up at low temperatures, while the U-vortex is not. Thus the vortex charge can be discussed in terms 
of the vortex bound states for all types of chiral pairings. Similarly, we can understand that the s-wave (zero winding) 
vortex is charged upJlj 

Finally let us discuss non-chiral pairings, such as p Xl p yir a.nd d x 2_ y 2-waves. These states are all gapless and it is 
reported that there is no bound states for the d x 2_ y 2-w&ve.cQ Since the system is anisotropic in this case, the energy 
eigenstates are not classified by the angular momentum. Therefore the above argument for the vortex charging is 
difficult to apply to this case. However, our result implies that the detailed quasiparticle structure in the vortex core 
is important for the vortex charging effect even in the gapless pairing cases. 

IV. SUMMARY AND DISCUSSIONS 

In conclusion, we have solved the problem of a single p x ±ip y -wave vortex self-consistently within the Bogoliubov-de 
Gennes theory. The full self-consistent calculation including both charge and current screenings was performed for the 
first time to investigate vortex problems. Due to the time reversal symmetry breaking, there are two types of vortices. 
We found a substantial vortex charge in the C-vortex case, while the vortex charge is suppressed for U-vortex. We 
conclude that the vortex charging effects are mainly determined by the local quasiparticle structure around the vortex 
core, reflecting the chiral p x ± ip y -wave paring. Especially the lowest vortex bound state is very important for the 
charging effect at low temperatures. 

In a real sample we can expect two types (p x ± ip y ) of domains. Therefore we expect both C and U-vortices in 
the presence of an external magnetic field. In domains, where the C-vortex is realized, an electronic field is induced 
due to the vortex charge. Therefore the charge of the C-vortex can be detected. If we use a field cooled sample, a 
single domain forming the U-vortex is realized all over the sample, since the U-vortex is energetically favorable Ed In 
this case it is difficult to find a signal from the U-vortex, since the charge of the U-vortex is much smaller than the 
C-vortex one. Thus we can distinguish the two types of C and U-vortices. 

How can we detect the vortex charge? The vortex charge can induce lattice distortion_and it scatters neutrons. 
It is reported that a polarized neutron scattering can be used to detect the vortex charge.EZI NQR is also one of the 
possible experiment which can detect the vortex charging effect, since the NQR detects the local electric field induced 
by the vortex change. Very recently, Kumagai et al. reported that the vortex charge is observed by the NQR in a 
high-T c material.E3 Thus the detection of the vortex charge is in progress now. It is very exciting if the vortex charge 
of the chiral superconductor is detected, since it can provide a very strong evidence of the rotating Cooper pair. 
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APPENDIX A: DERIVATION OF BOGOLIUBOV-DE GENNES EQUATION AND GAP EQUATION 

In this appendix we derive the Bogoliubov-dc Gennes equation and gap equation for the p x + ip y -w&ve supercon- 
ductor, which we use in the text. 

Let us start from the the following Bogoliubov-de Gennes equation: 



h (r)u n (r) + J dr'A(r,r')v n (r') = E n u n (r), (Ala) 
-h* (r)v n (r) - f dr'A*(r,r')u n {r') = E n v n (r). (Alb) 



Coordinates r and r' can be transformed into the center of mass and relative coordinates. 



R=^—, X = r r'. (A2) 



The p x + ipy-wave pairing is expressed as 



A(R,k)=A x (R)-^+iA y (R)-^-, (A3a) 



A(R,X) = f dke lkX A(R,k) 

= l (2^p j dk±-[A x (R)d x , + iA y (R)d y ,]e lk < r - r "> 



= -^-[A x (R)d x >+iA y (R)d y ,]5{r-r'). (A3b) 



Here A x and A y are the p x and p y components of the order parameter. Substituting Eq. (|A3|b) into Eq. (Al), we 
obtain 

h (r)u n (r) - ^-^A x (r)d x + iA y (r)d v + ^ [d x A x (r) + id y A v {r)] }v n (r) = E n u n (r), (A4a) 

-h *(r)v n (r) - ^-^A x (r)d x - iA y {r)d y + ~[d x A x (r) - id y A y {r)] }w n (r) = E n v n (r). (A4b) 
Next we derive gap equation. The order parameter is defined by 

A(r,r')=V(r,r') £ tanh u n (r')v n *(r) = V(r,r')D(r,r'), (A5) 

where E c is a cutoff energy. For the p-wave pairing, we assume the following p-wave interaction: 



V(r, r') = -j^- / dke — j. (A6) 



W (dke^ 2 

(2tt) 2 J k F 

Here V p > and Q are the p-wave attractive interaction and two-dimensional volume of the system, respectively. The 
order parameter is then expressed as 



A(R,k)= J dXe- lk x A(r ir ') 



= -V p n J dk'^j^ J dXe- tkX D(R +y R ~ Y )eMX - (A7a) 



We expand D{R + x -, R - ^) as 



X X, „, T dD(R,R') dD(R,R r ), X 

D{R+—,R-—) ~ R & R) + r'—>r • T . (A8) 

and obtain 

F 0<£„<B C ^ ' 

- lVpn ¥- { M-^ ] ^ tznh(^)u n (R)v:(R')\ R ^ R . (A9) 

F 0<B„<B c 

The first term in Eq. ( |A9| ) is zero due to the p-wave pairing. Since A(_R, fe) can be divided into two parts as in Eq. 
(|A3|a), the gap equation takes the following form: 

A x (r) = -V p fl^(d x - d x >) tanh ( ^ U„(rK*(r')|r'^r, (AlOa) 

F 0<E„<E c ^ ' 

\{r) = -V p Cl^-{d y - d yl ) £ t£Ulh (§) u n (r)v n *(r')\ r ,^ r . (AlOb) 

F 0<E„<E c ^ ' 

For a cylindrical system, it is convenient to introduce the following form for order parameters: 



A ± (r) = -[A x (r)±A»]. (All) 



Here A± is the order parameter for the p x ± p y -wave pairing. By using the A± representation, the Bogoliubov-de 
Gennes equation and the gap equations are expressed as 

h (r)u„(r) ~ ^-{A+(r)n+ + A_(r)n_ + -[a + A+(r) + Q_A_(r)] }u„(r) = E n u n (r), (A12a) 

-ftS(r)« n (r)-^{A + (r)n++A_(r)n_ + i[n + A + (r) + D_A_(r)]}V(r)=^w B (r), (A12b) 

A ± (r) = -^-^(D T -D' T ) £ taah (^W^K V)lw^-, (Al2c) 

U±=d x ±id y . (A12d) 
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FIG. 1. Self-consistent results at T — 0. (a) Order parameters scaled by Ao, which is the order parameter at T — in the 
bulk region. £o = vf/Aq is the coherence length. Set of parameters are taken as R — 8£o, &f£o = 16, Al = fo, Atf = 1/&F, 
E c — fi. Here Al and Atf are the London and Thomas-Fermi screening lengths, respectively, (b) Bound state energy spectrum 
for the C-vortex (circle) and U-vortex (X). is the angular momentum of the wave function u n . The bound states represented 
by Bci-ajid Bu are important in the vortex charging effect (see the text). Extended states lie continuously in a \E n \ > Ao 
region .Ej 
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FIG. 2. (a) Dimensionless scalar potential ao — eAo/Ao, electric field e r = —^od r ao, and charge density p = £o(l/r + d r )e r 
for the C- vortex at T = 0. (b) Spatial dependence of the electron density n in 1/£q unit at various temperatures. Charge 
density is given by p(r) — e[n(oo) — n(r)]. 




FIG. 3. Local electron density for the C- vortex in 1/£q unit, (a) Contribution from the bound state Be [Fig. 0(b)]. It is 
zero at T = 0. (b) Bound state contribution tib- 
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FIG. 4. Total electron density n and contribution from the extended states tie at T — 0. Solid lines represent the result 
with charge screening (A T p = 1/&f), while dashed lines are the result without charge screening (A T p — » oo). 
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FIG. 5. Schematic picture of bound state energy spectrum for various winding pairings. Numbers represent the angular 
momenta of u n and v n (l u and l v ). 
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